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ABSTRACT 

The long-held notion that the highest-energy cosmic rays are of distant 
extragalactic origin is challenged by observations that events above ~ 10 20 
eV do not exhibit the expected high-energy cutoff from photopion production 
off the cosmic microwave background. We suggest that these unexpected 
ultra-high-energy events are due to iron nuclei accelerated from young strongly 
magnetized neutron stars through relativistic MHD winds. We find that neutron 
stars whose initial spin periods are shorter than ~ 10 ms and surface magnetic 
fields are in the 10 12 — 10 14 G range can accelerate iron cosmic rays to greater 
than ~ 10 20 eV. These ions can pass through the remnant of the supernova 
explosion that produced the neutron star without suffering significant spallation 
reactions or energy loss. For plausible models of the Galactic magnetic field, the 
trajectories of the iron ions curve sufficiently to be consistent with the observed, 
largely isotropic arrival directions of the highest energy events. 

Subject headings: acceleration of particles, magnetic fields, MHD, plasmas 



The detection of cosmic rays with energies above 10 20 eV has triggered considerable 
interest on the origin and nature of these particles. Hundreds of events with energies 
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above 10 19 eV and about 20 events above 10 20 eV have now been observed by a number of 
experiments such as HiRes (Kieda et al. 1999), AGASA (Takeda et al. 1998, 1999), Fly's 
Eye (Bird et al. 1995, 1993, 1994), Haverah Park (Lawrence, Reid & Watson 1991), Yakutsk 
(1996), and Volcano Ranch (1963). Most unexpected is the large flux of events observed 
above 5 x 10 19 eV (Takeda et al. 1998) with no sign of the Greisen-Zatsepin-Kuzmin (GZK) 
cutoff (Greisen 1966; Zatsepin & Kuzmin 1966). The cutoff should be present if these 
ultra-high energy particles are protons produced by sources distributed homogeneously 
throughout the universe. Cosmic ray protons of energy above 5 x 10 19 eV lose their energy 
to photopion production off the cosmic microwave background and cannot originate further 
than about 50 Mpc away from us. Alternatively, if ultra-high-energy cosmic rays (UHECRs) 
are protons from sources closer than 50 Mpc, the arrival direction of the events should point 
toward their source. The present data shows a mostly isotropic distribution and no sign of 
the local distribution of galaxies or of the Galactic disk above 10 19 eV (Takeda et al. 1999). 
In sum, the origin of these particles with energies tens of millions of times greater than any 
produced in terrestrial particle accelerators, remains a mystery. 

In addition to the difficulties with locating plausible sources of UHECRs in our nearby 
universe, there are great difficulties with finding plausible accelerators for such extremely 
energetic particles. Acceleration of cosmic rays in astrophysical plasmas occurs when the 
energy of large-scale macroscopic motion, such as shocks and turbulent flows, is transferred 
to individual particles. The maximum possible energy, E max , is estimated by requiring that 
the gyro-radius of the particle be contained in the acceleration region (Hillas 1994) and 
that the acceleration time be smaller than the time for energy losses. The former condition 
relates E max to the strength of the magnetic field, B, and the size of the acceleration 
region, L, such that -E max ~ ZeB L, where Ze is the charge of the particle. For instance, 
for -Emax ~ 10 20 eV and Z ~ 1, the known astrophysical sources with reasonable BL 
products are neutron stars (B ~ 10 12 G and L ~ 10 km), active galactic nuclei (B ~ 10 4 
G and L ~ 10 AU), radio galaxies {B ~ 10~ 5 G and L ~ 10 kpc), and clusters of galaxies 
(B ~ 10~ 6 G and L ~ 100 kpc) (Hillas 1984; Berezinsky et al. 1990). However, energy 
losses usually prevent acceleration to E mSLX , and no effective mechanism for UHECR 
acceleration has been shown for any of these objects (Blandford 1999; Bhattacharjee & Sigl 
1998; Venkatesan, Miller & Olinto 1997). Here we show that the early evolution of young 
magnetized neutron stars in our Galaxy may be responsible for the flux of cosmic rays 
beyond the GZK cutoff. A preliminary study of this idea can be found in (Olinto, Epstein 
& Blasi 1999). 

Neutron stars have been previously suggested as possible sources of UHECRs, starting 
with an early attempt by Gunn & Ostriker (1969) to the more recent proposal of Bell 
(1992). Thus far, these attempts have failed at either reaching the highest energies or 



- 3- 



reproducing the spectrum or the apparent isotropy of the arrival directions of UHECRs. 
In the following we describe our alternative. We propose that young neutron stars may 
accelerate heavy nuclei to the highest observed energies by transferring their rotational 
energy to particle kinetic energy via a relativistic MHD wind. 

Some neutron stars may begin their life rotating rapidly (Q ~ 3000 rad s _1 ) and with 
large surface magnetic fields (B$ ^ 10 13 G). The dipole component of the field decreases 
as the cube of the distance from the star's surface B(r) = BsiRs/r) 3 , where the radius of 
the star is R,s — 10 6 cm. As the distance from the star increases, the dipole field structure 
cannot be causally maintained, and beyond the light cylinder radius, R\ c = c/Q, the field 
is mostly azimuthal, with field lines spiraling outwards (Michel 1991). For young, rapidly 
rotating neutron stars, the light cylinder is about ten times the stellar radius, R[ c = 10 7 Q 3l ^ 
cm, where f2 3fe = fi/3000 rad s -1 . 

The surface of young neutron stars is composed of iron peak elements formed during 
the supernova event. Iron ions can be stripped off the hot surface of a young neutron 
star due to strong electric fields and be present throughout much of the magnetosphere 
(Ruderman & Sutherland 1975, Arons & Scharlemann 1979). Inside the light cylinder, 
the magnetosphere corotates with the star and the iron density has the Goldreich- Julian 
value: naj{r) = B (r)Q / (4tt Z ec) , where c is the speed of light (Goldreich & Julian 1969). 
In this estimate, and what follows, we do not include the trigonometric factors related to 
the relative orientation of the magnetic and rotational axes. 

The exact fate of the plasma outside the light cylinder is still a subject of debate 
(Gallant k Arons, 1994; Begelman & Li, 1994; Chiueh, Li, & Begelman, 1998; Melatos 
& Melrose, 1996). Observations of the Crab Nebula indicate that most of the rotational 
energy emitted by the Crab pulsar is converted into the kinetic energy of particles in a 
relativistic wind (Kennel & Coroniti 1994; Begelman 1998; Emmering & Chevalier 1987). 
This conversion may be due to properties of the MHD flow, related to magnetic reconnection 
(Coroniti 1990), or a more gradual end of the MHD limit (Melatos & Melrose, 1996). Some 
analytical and numerical studies show the development of kinetically dominated relativistic 
winds (see e.g., Begelman and Li, 1994), but at present the theoretical understanding of the 
wind dynamics is far from complete. 

The basic idea of accelerating plasmas by the Poynting flux was proposed by Weber 
and Davis (1967) (then called magnetic slingshot). Later, Michel(1969) showed that for 
a perfectly spherical flow the complete conversion of the magnetic energy into kinetic 
energy of the flow could not be achieved. However, Begelman and Li (1994) reconsidered 
the problem and showed that even small deviations from a spherical flow could imply an 
efficient conversion of the magnetic energy into kinetic energy of the wind through the 
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so-called magnetic nozzle effect, provided the magnetic field lines have the right geometry. 

In the present study we assume that, at least for some neutron stars most of the 
magnetic energy in the wind zone is converted into the flow kinetic energy of the particles 
in the wind and that the rest mass density of the regions of the wind containing iron ions 
are not dominated by electron-positron pairs; that is, the electron-positron density is less 
than ~ 10 5 times that of the iron ions. With these assumptions, the magnetic field in the 
wind zone decreases as B(r) <^ Bi c Ri c /r. For surface fields of B$ = 10 13 -Bi3 G, the field at 
the light cylinder is Bi c = 10 10 B\ 3 Q 3k G. The maximum energy of particles that can be 
contained in the wind near the light cylinder is 

E max = ^1^1 ~ 8 x 10 20 Z 26 B 13 nj k eV , (1) 
c 

where Z 2 & = Z/26. In the rest frame of the wind, the plasma is relatively cold while in the 
star's rest frame the plasma moves with Lorentz factors ~ 10 9 — 10 10 . 

The typical energy of the accelerated cosmic rays, E cr , can be estimated by considering 
the magnetic energy per ion at the light cylinder E cr ~ BfJ^imcj- At the light cylinder 
n GJ = 1.7 x 10 11 B 13 ttj k /Z cm~ 3 which gives 

E cr ~ 4 x 10 20 Z 26 B 13 Q 2 3k eV , (2) 

similar to E max above (Gallant & Arons 1994; Begelman 1994). 

The spectrum of accelerated UHECRs is determined by the evolution of the rotational 
frequency: As the star spins down, the energy of the cosmic ray particles ejected with the 
wind decreases. The total fluence of UHECRs between energy E and E + dE is 

where the particle luminosity is 

N = £ n G j nRlc = 6 x i 3 ^^!l B -i (4 ) 

and £ < 1 is the efficiency for accelerating particles at the light cylinder. The rotation 
speed decreases due to electromagnetic and gravitational radiation (Lindblom, Owen & 
Morsink 1998; Andersson, Kokkotas & Schutz 1999). For B s ^ 10 13 G, r-mode gravitational 
radiation is likely suppressed (Rezzolla, Lamb & Shapiro 1999) and the spin down may be 
dominated by magnetic dipole radiation given by: 

It* = (5) 
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For a moment of inertia I = 10 45 g cm 2 , the time derivative of the spin frequency is, 
Cl= 1.7 x lO-^fijjfcS- 1 , and Eq. (|) gives 

Substituting in Eq. (||), the particle spectrum from each neutron star is 

N(E) = e^^GeV- 1 , (7) 

where E = 10 20 E 20 eV. 

Neutron stars are produced in our Galaxy at a rate 1/r, where r = 100 r 2 yr, and a 
fraction e of them have the required magnetic fields, initial spin rates and magnetic field 
geometry to allow efficient conversion of magnetic energy into kinetic energy of the flow. 
As discussed below, UHE iron nuclei scatter and diffuse in the Galactic magnetic field. 
Taking the confining volume for these particles to be V c and the lifetime for confinement to 
be t c , the UHECR density is n(E) = eN(E)t c /rV c , and the flux at the surface of the Earth 
is F(E) = n(E)c/4. For a characteristic confinement dimension of R = 10 Ri kpc we can 
write V c = AttR 3 /3 and t c = QR/c, where Q > 1 is a measure of the how well the UHECR 
are trapped. The predicted UHECR flux at the Earth is 

F(E) = 10' 24 r p J €Q F ? GeV-WV 1 . (8) 

By comparing with observations, we can estimate the required efficiency 
factor, £e. The AG AS A experiment finds that the flux at 10 20 eV at Earth is 
F{E) = 4 x 10^ 30 GeV^ 1 cm" 2 s~ 1 . Equating this flux with the estimate of Eq. (D, 
we find that the efficiency factor only needs to be £e ^ 4 x 10~ 6 Q -1 - The smallness of 
the required efficiency suggests that young, Galactic neutron stars can be the source of 
UHECRs even if only a small fraction of stars are born with very rapid spin frequencies and 
high magnetic fields. 

The observed energy spectrum of cosmic rays below the expected GZK cutoff (i.e., 
between ~ 10 8 eV and ^ 10 19 eV) has a steep energy dependence N(E) oc E' 1 , with 7 m 2.7 
for E £ 10 15 eV and 7 w 3.1 for 10 15 £ E(eV) £ 10 19 (Gaisser 1990). The events with 
energy above 10 19 ' 5 eV, however, show a much flatter spectrum with 1 ^ 7 <^ 2; the drastic 
change in slope suggests the emergence of a new component of cosmic rays at ultra-high 
energies. The predicted spectrum of Eq. (|8|) is very flat, 7 = 1, which agrees with the 
lower end of the plausible range of 7 observed at ultra-high energies. Propagation effects 
can produce an energy dependence of the confinement parameter Q and, correspondingly, a 
steepening of the spectrum toward the middle of the observed range 1 ^ 7 <; 2. 
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Even though a young neutron star is usually surrounded by the remnant of the 
presupernova star, the accelerated particles can escape the supernova remnant without 
significant degradation as the envelope expands. A requirement for relativistic winds to 
supply UHECRs is that the column density of the envelope becomes transparent to UHE 
iron nuclei before the spinning rate of the neutron star decreases to the level where the star 
is unable to emit particles of the necessary energy. 

To estimate the evolution of the column density of the envelope, consider a supernova 
that imparts Esn = 10 51 £si erg to the stellar envelope of mass M env = 10 Mi M . The 
envelope then disperses with a velocity v e ~ {2E SN / M env ) 1 ^ 2 = 3 x 10 8 (S^/Mi) 1 ^ 2 cm s -1 . 
The column density of the envelope surrounding the neutron star is S ~ M env jAv:R 2 e ^ 
where R e ff = Ro + v e t, where Ro is the characteristic radius of the presupernova star, 
Rq <; 10 14 cm. We now have 

M env = 16 cm - 2 (9) 

where t is in seconds, and t e = Ro/v e ^3 x 10 5 (Mi/£ 5 i) 1 / 2 s. The condition for iron nuclei 
to traverse the supernova envelope without significant losses is that £ <^ 100 g cm -2 . This 
"transparency" occurs at times t > t tr = 1.3 x 10 7 M 1( ? 51 1/2 s > t e . 

As the envelope is being ejected, the neutron star spin is slowing due to the magnetic 
dipole radiation, Eq. (|5|), so that 

where 3000 rad s -1 is the initial spin rate and tg — t/10 s s. The cosmic ray energy thus 
evolves according to 

E cr (t) = 4 x 10 2 °eV n Z lffl" . (11) 

[1 + *8-Dl3"i3A:J 

The condition that a young neutron star could produce the UHECRs is that E cr exceeds 
the needed energy when the envelope becomes transparent; i.e., E cr (t tr ) > 10 20 E 2 o eV. This 
translates into the following condition: 



3000 s- 1 

> — ~r ^ • (12) 

4Z 26 E 20 1 -0.13M 1 5 13 ^ 5 " 1 1/21 
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From this equation we obtain the allowed regions in the Bs~Qi plane shown in Figure 1 for 
E 2 q = 1 and 3 and M env = 5 and 50 M Q . 

For the parameters within the allowed region, the acceleration and survival of UHE iron 
nuclei is not significantly affected by the ambient photon radiation. The most important 
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source of radiation in the wind region is the thermal emission from the star's surface. The 
low energy non-thermal radiation from the neutron star is not significant unless it is > 10 4 
that of the Crab pulsar. In the time needed for the envelope to become transparent, the 
surface cools to ~ 3 x 10 6 K (Tsuruta 1998). For these temperatures, photodissociation 
(see, e.g., Protheroe, Bednarek & Luo 1998) and Compton drag have minor effects on the 
energy and composition of the accelerating iron nuclei. Furthermore, synchrotron losses are 
unimportant because the plasma is essentially cold in the rest frame of the accelerating 
plasma. 

The relativistic MHD wind from a rapidly spinning neutron star may impart more 
energy to the supernova remnant than the initial explosion. For initial spin rates ^ 1000 rad 
s _1 , the rotational energy is ^ 10 51 erg, comparable to the kinetic energy of most supernova 
remnants. More rapidly spinning neutron stars may generate highly-energetic supernova 
events, possibly similar to SN 1998bw (Kulkarni 1998). In these cases, the right boundary 
of the allowed region in Figure 1 should be enlarged because the remnant expands more 
rapidly than assumed above. 

The iron ejected with energies ~ 10 20 eV will reach Earth after being deflected by the 
Galactic and halo magnetic fields (Zirakashvili, Pochepkin, Ptuskin & Rogovaya 1998). The 
gyroradius of these UHECRs in the Galactic field of strength B gai is 



which is considerably less than the typical distance to a young neutron star (~ 8 kpc). 
Therefore, ultra-high energy iron arriving at the Earth would not point at the source. 
A Galactic iron source is consistent with an approximately isotropic arrival direction 
distribution as observed by AGASA for UHECRs (Zirakashvili et al. 1998). In support of 
this interpretation, we note that the cosmic ray component at 10 18 eV is nearly isotropic 
with only a slight correlation with the Galactic disk and spiral arms (Hayashida et al. 
1999). If these cosmic rays are protons of Galactic origin, their isotropy is indicative of the 
diffusive effect of the Galactic and halo magnetic fields. Since the iron arrival distribution 
at 10 20 eV probes similar trajectories to protons at a few times 10 18 eV we expect the iron 
to show a nearly isotropic distribution with a slight correlation with the Galactic center 
and disk. This correlation should become apparent if the number of observed events grows 
by orders of magnitude or if events with energies higher than the present highest energies 
events are detected. Although some indication of a correlation with the Galactic center for 
events above 10 20 eV has been recently reported (Stanev & Hillas 1999), the small number 
of observed events limits the significance of this finding. 
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In conclusion, we propose that ultra-high-energy cosmic ray events originate from 
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iron nuclei accelerated by young, strongly magnetic, Galactic neutron stars. Iron from the 
surface of newborn neutron stars are accelerated to ultra-high energies by a relativistic 
MHD wind. Neutron stars whose initial spin periods are shorter than ~ 4(£?s'/10 13 G) 1 / 2 ms 
can accelerate iron nuclei to greater than 10 20 eV. These ions can pass through the radiation 
field near the neutron star and the remnant of the supernova explosion that produced the 
neutron star without suffering significant deceleration or spallation reactions. 

The best test of this proposal is a unambiguous composition (mass/charge) 
determination and a correlation of arrival directions for events with energies above 10 20 eV 
with the Galactic center and disk. Both aspects will be well tested by future experiments 
such as the Auger Project (Cronin 1999) and OWL-Airwatch (Ormes et al., 1997). In 
addition, our model will be severely constrained if the indication of a small scale clustering 
among UHECR events (Uchihori et al. 1999) is confirmed by future experiments to be due 
to an isotropically distributed set of discrete sources. 
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Figure caption 

Fig. 1 

Parameter space for which acceleration and escape of the accelerated particles through 
the ejecta are allowed. The solid lines refer to particle energy E cr = 10 20 eV and dashed 
lines to E cr = 3 x 10 20 eV. The curves are plotted for two values of the envelope mass, 
M env = 50 M so i ar and M env = 5 M so i ar , as indicated. The horizontal line at spin period 
~ 0.3 ms indicates the minimum period (maximum angular speed) allowed for neutron stars 
(Haensel, Lasota & Zdunik 1999) . 



